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FOLD-SADDLE BIFURCATION IN NON-SMOOTH 
VECTOR FIELDS ON THE PLANE. 

CLAUDIO A. BUZZI\ TIAGO DE CARVALHO^ AND 
MARCO A. TEIXEIRA^ 



Abstract. This paper presents results concerning bifurcations of 2D 
piecewise— smooth dynamical systems governed by vector fields. Generic 
three— parameter families of a class of Non— Smooth Vector Fields are 
studied and its bifurcation diagrams are exhibited. Our main results 
^^ ' describe the unfolding of the so called Fold — Saddle singularity. 

q 

^ ' 1. Introduction 

cd . 

The general purpose of this article is to present some aspects of the geo- 
metric and qualitative theory of a class of planar non— smooth systems. Our 
main concern is to discuss the behavior of such systems around typical sin- 
gularities that appear generically in three— parameter families. We mention 
00 ' that certain phenomena in control systems, impact in mechanical systems 

lO . and nonlinear oscillations are the main sources of motivation of our study 

^^ I concerning the dynamics of those systems that emerge from differential equa- 

tions with discontinuous right— hand sides. 

The codimension zero and codimension one singularities were discussed 
(•^ I in [4] and [5] respectively. In [3] codimension two singularities were studied. 

The specific topic addressed in this paper is the complete characterization 

of the Fold— Saddle bifurcation diagram. Those papers give the necessary 

(s^ ■ basis for the development of our approach. 

^; 

Cd I Let i^ C M^ be a compact set and S C X given by S = / ^(0), where 

/ : iiT — >• M is a smooth function having G M as a regular value (i.e. 
Vf{p) ^ 0, for any p G /~^(0)) such that dK nT, = 9 or dK (h S. Clearly 
T, is the separating boundary of the regions T,^ = {q ^ ^IfiQ) ^ 0} s-nd 
S_ = {q G K\f{q) < 0}. We can assume that S is represented, locally 
around a point q = (x, y), by the function f{x, y) = y. 

Designate by x^ the space of C' vector fields on K endowed with the 
C^— topology with r > 1 or r = cx), large enough for our purposes. Call 
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n'' = n''{K, f) the space of vector fields Z : K \ S ^ M^ such that 



Z{x,y) 



X{x,y), for (x,y) S S+, 
Y{x,y), for (x,y)GS_, 



where X = {fi,gi), Y = (/2,5'2) are in x^- We write Z = {X,Y), which 
we will accept to be multivalued in points of S. The trajectories of Z are 
solutions of q = Z{q), which has, in general, discontinuous right— hand side. 
The basic results of differential equations, in this context, were stated by 
Filippov in [2]. Related theories can be found in [HEIIH]. 

In what follows we will use the notation 

X.fip) = {Vf{p),X{p)) and Y.fip) = {Vf{p),Y{p)) . 



1.1. Setting the problem. Let Xq be a smooth vector field defined in S+. 
We say that a point po £ ^ is a T,—fold point of Xq if Xo./(po) = but 
Xq./(po) 7^ 0. Moreover, po G 5] is a visible (respectively invisible) S— fold 
point of Xo if Xo./(po) = and X^.f{po) > (resp. X^.f{po) < 0). In this 
universe, T^", a S— fold point has codimension zero. Since f{x,y) = y we 
derive the following generic normal forms Xq{x, y) = (ai, /3ix) with ai = ±1 
and /3i = ±1. 

Let Yq be a smooth vector field defined in S_. Assume that Yq has a 
hyperbolic saddle point Syq on S and that the eigenspaces of DYq{Syo) 
are transverse to S at Syq- In this universe, '^•^, a saddle point Syq has 
codimension one. Since f{x,y) = y we derive the following generic normal 
forms YQ{x,y) = {a2y,a2x) with 02 = ±1 and its generic unfolding Yjs = 
(02(2/ + /3), a2x) where /3 € M. Let U he a small neighborhood of Yq in F^,". 
Then: 

(a) There exists a smooth function L : U ^ M, such that DLyq is surjective. 

(b) The correspondence Y -^ Sy is smooth, where Sy is a saddle point of 
Y. 

(c) If L(y) >0 thenar € S_. 

(d) If L{Y) = then Sy € S. 

(e) IfL(y) <0then5y G S+. 

In this paper we are concerned with the bifurcation diagram of systems 
Zq = {Xq,Yq) in Q^' such that po = Syq G S. This singularity will be called 
Fold — Saddle singularity (see Figures [1] and [2]) . 

We depart from Zq,Zq G il.^ written in the following forms: 



(1) 



ZA 



X^ 



Yq 



1 

-x 

-y 

-X 



if y > 0, 

if y < 0, and 
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Figure 1. (Invisible) 
Fold- Saddle Singularity. 



Figure 2. (Visible) 
Fold-Saddle Singularity. 



(2) 



Zv 




^0 



^n 



if y > 0, 
if y < 0. 



Note that X'q presents an invisible S— fold point on its phase portrait and 
Xq presents a visible one. Following the techniques developed in [^, we are 
able to prove that there exists a smooth mapping F-r : Q^',Zq — ?• M'^, where 
T = i or V such that: 

1- {DFr)z^ is surjective (So Mr = {Fr)^^{0) is locally, around Zq, an 
imbedded differentiable manifold). 

2- Each Z £ Ut, with Fr{Z) = and Ut a small neighborhood of Zg in 
O^ is C"— equivalent to Zq. 

The main question is to exhibit the bifurcation diagram of Zq. So, we have 
to consider generic imbeddings cJt- : M'^, — > fi'', Zq (3— parameter families). 
They are transversal imbeddings to Mr at Zq. 

Consider Zq = {Xq.,Yq) G Ur. Roughly speaking, we derive that: 
I- There is a canonical imbedding Fq : 'S? ,Q ^ x^ ■, ^q such that Fq{X, (3) = 
ZJ^ a expressed by: 



(3) 



^X,I3 



XI 



Yr 



-(2/ + /3) " 



if y > 0, 
if y < 0, 



where A, /3 G (—1, 1), a\{i) = —1 and a\{v) = 1. Moreover, its bifurcation 
diagram of Zj^ a is exhibited (see Figures [18] and [28l) . We observe that there 
are some typical topological types nearby Zq that do not appear in the 
bifurcation diagram of ZJ o. For example, when t = i the configurations in 
Figures [3] and [4] are excluded and when r = v the configuration in Figure [5] 
also is excluded. 
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Figure 3. 



Figure 4. 



Figure 5. 



II- We add an auxiliary parameter // in the following way: 



(4) 



Z 



A,M,/3 



x^ 



1 



ai(r)(a; — A) 



Y, 



M,/3 



I 



X + 
>-2) 



(/.,-2) 



(y + /3) 



^ + f(y + /3) 



if y > 0, 
if y < 0, 



where A, /3 G (—1, 1), Qi(i) = —1, a\{v) = 1 and fi € (— eoi^o) with the real 
number eo > being sufficiently small. By means of this late unfolding its 
bifurcation diagram cover all topological types near ^o co- 
in this universe, F^o ^ a Fold— Saddle singularity has codimension three. 
Since f{x,y) = y we derive the generic normal forms ^a,/xo,/3 with /io = 
±eo/2 and its generic unfolding Z^^^^p given by (j4]). Therefore, there is a 
codimension three bifurcation (global) branch terminating at Zq. In fact, 
note that we can obtain Equation ([1]) (respectively dJ])) from Equation @ 
taking t = i (respectively r = f),A = 0, /i = and /3 = . 

Of course, we can take another generic normal form of one or both vector 
fields Xq and Yq- Iii this paper we consider just the cases described in 
Equations ([TJ and ([2]). For the other cases a similar approach can be done. 
It is worth mentioning that we detect branches of ^^ canard cycles^^ in the 
bifurcation diagram of Z^^^j^. Recall that, a canard cycle is a closed path 
composed by pieces of orbits of X, Y and Z^ (see Figures [71 [8] and [9]). In 
Section 2 a precise definition will be given. 

Example 1. Equations ([T]) and ([2]) appear in problems related to Control 
Theory, more specifically, in Relay Systems. In fact, consider the function 
(/3 : M — >• M given by 

-1, fory>0, 

-y, for y <0, 



v>{y) 



and u{y) = —ip{y)sign{y). So ([T]) and ([2]) are represented by Z^^x^y) = 
{u{y),a{T)x) where t = i or v, a{i) = —1 and a{v) = sign{y). 

1.2. Statement of the Main Results. Our results are now stated. The- 
orems 1, 2 and 3 are intermediate steps towards Theorem A and Theorems 
4, 5 and 6 are intermediate steps towards Theorem B. 
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Theorem 1. Take t = i in Equation ([3]) or equivalently, take t = i and fi = 
in Equation (j4]). The {X, f3)— plane contains essentially 17 distinct typical 
configurations representing 5 distinct topological behaviors on its bifurcation 
diagram (see Fiaure[TE\). 



It is easy to see that the cases covered by Theorem 1 do not represent the 
full unfolding of the (Invisible) Fold— Saddle singularity. Because of this, 
the next two theorems are necessary. Each one of them describes a distinct 
generic co dimension two singularity. 

Theorem 2. Take t = i and < fi < Eq in Equation (j4]). The {X, /3)— plane 
contains essentially 19 distinct typical configurations representing 7 distinct 
topological behaviors on its bifurcation diagram (see Fiaure[2U\). 



Theorem 3. Take t = i and —Eq < /i < m Equation @. The 
(A, (3)— plane contains essentially 19 distinct typical configurations represent- 
ing 7 distinct topological behaviors on its bifurcation diagram (see FiaureW^). 



Theorem 4. Take t = v in Equation ([3]) or equivalently, take t = v and 
fj, = in Equation ^. The {X, I3)~plane contains essentially 13 distinct 
typical configurations representing 7 distinct topological behaviors on its bi- 
furcation diagram (see Figure [26 



The cases covered by Theorem 4 do not represent the full unfolding of the 
(Visible) Fold— Saddle singularity. Because of this, the next two theorems 
are necessary. Each one of them describes a distinct generic codimension 
two singularity. 

Theorem 5. Take t = v and < /i < Eq in Equation Q. The (A, P)— plane 
contains essentially 13 distinct typical configurations representing 7 distinct 
topological behaviors on its bifurcation diagram (see Figure\ 



Theorem 6. Take t = v and —Eq < fi < in Equation ^. The 
{X, P)— plane contains essentially 13 distinct typical configurations represent- 
ing 7 distinct topological behaviors on its bifurcation diagram (see Fiaure[28\). 

Finally, we are able to state the main results of the paper. 

Theorem A. Equation Q with t = i generically unfolds the (Invisible) 
Fold— Saddle singularity. Moreover, its bifurcation diagram exhibits 55 dis- 
tinct typical configurations representing 11 distinct topological behavior (see 
FigureWB(l. 
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Theorem B. Equation ^ with t = v generically unfolds the (Visible) 
Fold— Saddle singularity. Moreover, its bifurcation diagram exhibits 39 dis- 
tinct typical configurations representing 21 distinct topological behavior (see 
Figure\3C 



The paper is organized as fohows: in Section 2 we give the basic theory 
about Non— Smooth Vector Fields on the Plane, in Section 3 we prove The- 
orem 1, in Section 4 we prove Theorem 2, in Section 5 we prove Theorem 
3, in Section 6 we prove Theorem A and present the Bifurcation Diagram 
of Z)^ n fl; in Section 7 we prove Theorem 4, in Section 8 we prove Theorem 
5, in Section 9 we prove Theorem 6 and in Section 10 we prove Theorem B 
and present the Bifurcation Diagram of Z^^ o. 

2. Preliminaries 

We distinguish the following regions on the discontinuity set S : 
(i) Si C S is the sewing region if {X.f){Y.f) > on Si . 
(ii) S2 C S is the escaping region if {X.f) > and (Y.f) < on S2. 

(iii) S3 C S is the sliding region if (X.f) < and (Y.f) > on S3. 

Consider Z £ Q^ . The sliding vector field associated to Z is the vector 
field Z^ tangent to S3 and defined at g G S3 by Z^{q) = m — q with m being 
the point where the segment joining q + X[q) and q + Y{q) is tangent to S3 
(see Figure [6]). It is clear that if (7 G S3 then (7 E S2 for —Z and then we 
can define the escaping vector field on S2 associated to Z hj Z'^ = —[—Zy. 
In what follows we use the notation Z^ for both cases. 




\q + X{q) 

Figure 6. Fillipov's convention. 

We say that (7 € S is a Tj— regular point if 
(i) {X.f{q)){Y.f{q)) > or 

(ii) {X.f{q)){Y.f{q)) < and Z^{q) ^ (that is g G S2 U S3 and it is 
not a singular point of Z^). 

The points of S which are not S— regular are called Ti — singular. We dis- 
tinguish two subsets in the set of S— singular points: S* and S^. Any g G S^ 



FOLD-SADDLE BIFURCATION IN NON-SMOOTH SYSTEMS 7 

is called a pseudo equilibrium of Z and it is characterized by Z (q) = 0. Any 
g G S* is called a tangential singularity and is characterized by Z (q) ^ 
and X.f{q)Y.f{q) = (g is a contact point of Z^). 

A pseudo equilibrium q € S^ is a T,— saddle provided one of the following 
condition is satisfied: (i) g G S2 and q is an attr actor for Z^ or (ii) g € S3 
and g is a repeller for Z^. A pseudo equilibrium g G S^ is a Ti—repeller 
(resp. T^—attractor) provided g^ G S2 (resp. q G S3) and q is a repeller 
(resp. attractor) equilibrium point for Z^. 

Definition 1. Consider Z ^Vf . 

(1) A curve T is a canard cycle ifV is closed and 

• r contains arcs of at least two of the vector fields -^|s+; ^|s_ 
and Z^ or is composed by a single arc of Z^; 

• the transition between arcs of X and arcs ofY happens in sewing 
points; 

• the transition between arcs of X (or Y) and arcs of Z hap- 
pens through T,—fold points or regular points in the escape or 
sliding arc, respecting the orientation. Moreover ifT^T, then 
there exists at least one visible T,—fold point on each connected 
component o/T fl S. 

(2) Let T be a canard cycle of Z. We say that 

• T is a canard cycle of kind I if T meets S just in sewing 
points; 

• T is a canard cycle of kind II seV = S; 

• T is a canard cycle of kind III if T contains at least one 
visible Ti—fold point of Z . 

In Figures [^ and [2| arise canard cycles of kind I, II and III 
respectively. 

(3) Let T be a canard cycle. We say that T is hyperbolic if 

• T is of kind I and ri'{p) ^ 1, where rj is the first return map 
defined on a segment T with p G T iti 7; 

• T is of kind II; 

• T is of kind III and or F n S C Si U S2 or F n S C Si U S3. 

Remark 1. The expression "canard" is used here because these orbits are 
limit periodic sets of singular perturbation problems (see [l]j. 

Definition 2. Consider Z G 0,^ . A point q (^ T, is a Ti— center if there is 
a neighborhood U of q such that an one parameter family of canard cycles 
encircles q and foliates U. 

Definition 3. Consider Z G^ . A closed path A is a T,— graph if it is a 

union of equilibria, pseudo equilibria, tangential singularities of Z and arcs 
of Z joining these points in such a way that A n S / 0. Like for canard 
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r<t 





Figure 


Figure 


Figure 


7. Canard 


8. Canard 


9. Canard 


cycle of 


cycle of 


cycle of kind 


kind I. 


kind II. 


III. 



cycles, we say that A is a T,— graph of kind I if A D'S C Si, A is a 
Ti— graph of kind II if An T, = A and A is a T,— graph of kind III if 

Ans c S2US3. 



In what follows, in order to simplify the calculations, we take /x = a + 1 
in ([3]) and obtain the following expression 



(5) 



'X,a,fS 



Xx 



^a,/3 



1 

ai(r)(x — A) 






-{y + P) 



'-X + 



if y > 0, 
if y < 0, 



\ 2 -^ ^ 2 

where A, /3 € (—1, 1), a € (—1 — eo, —1 + £0)1 t = i or v, ai{i) = —1 and 
ai(v) = 1. When it does not produce confusion, in order to simplify the 
notation we use Z = {X,Y) or Zx^a,/3 = {X,Y) instead -^Jq, « = (Xa)^«,^)- 

Given Z = {X,Y), we describe some properties of both X = Xx and 

The real number A measures how the S— fold point d = (A, 0) of X is 
translated away from the origin. More specifically, if A < then d is trans- 
lated to the left hand side and if A > then d is translated to the right hand 
side. 



Some calculations show that the curve Y.f = is given by y 



(i-Q) , 



/3. So the points of this curve are equidistant from the separatrices when 
a = — 1. It become closer to the stable separatrix of the saddle point 
S = Sa,(s when a G (—1, — l+eo). It become closer to the unstable separatrix 
of S when a S (— 1 — £9,— 1). Moreover, the smooth vector field Y has 
distinct types of contact with S according with the particular deformation 
considered. In this way, we have to consider the following behaviors: 

• Y~ : In this case /3 < 0. So S* is translated to the y— direction with 
y > (and S is not visible for Z). It has a visible S— fold point 

e = ea,(3 = (|Sy/^'0) (see Figure [TO]). 

• Y° : In this case /3 = 0. So 5 is not translated (see Figured]). 
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Y^ : In this case /? > 0. So 5 is translated to the y— direction 
with y < 0. It has an invisible S— fold point i = (ii,i2) = ^o,/3 = 






/3, ) . Moreover, we distinguish two points: h = hf 



-/3,0) 



which is the intersection between the unstable separatrix with S 
and j = jp = (/3, 0) which is the intersection between the stable 
separatrix with S (see Figure [TT]) . 

In Figure [TT] we distinguish the arcs cri of Y joining the saddle point S of 
Y to h and (T2 of Y joining j to the saddle point S oiY . 





Y.f:= 



Figure 10. Case Y~ 




Figure 11. Case Y+. 



3. Proof of Theorem 1 

In (a, 6) C E2 U S3, consider the point c = (ci,C2), the vectors X{c) = 
{di,d2) and Y{c) = (61,62) (as illustrated in Figure [T2|) . The straight seg- 
ment passing through c + X{c) and c + Y{c) meets S in a point p{c). We 
define the C— map 

p: {a,b) — y S 
z I — > p{z). 
We can choose local coordinates such that S is the x— axis; so c = (ci,0) 
and p{c) € M X {0} can be identified with points in M. According with this 
identification, the direction function on T, is defined by 

H : {a,b) — > R 

z I — > p{z) — z. 




c + X(c) 
Figure 12. Direction function. 



We obtain that H is a C*^'— map and 

• if H{c) < then the orientation of Z^ in a small neighborhood of c 
is from 6 to a; 
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(6) 



• if H{c) = then c G S^; 

• if H{c) > then the orientation of Z in a smah neighborhood of c 
is from a to b. 

Simple calculations show that p{ci) = '^^^ ^ '^l _/ &iicl consequently, 

62^1 - d^ei 



H{ci) 



62 - d2 

We now in position to prove Theorem 1. 

Proof of Theorem 1. In Cases li, 2i and 3i we assume that Y presents the 
behavior Y~ . In Cases 4i, 5i and 6i we assume that Y presents the behavior 
Y^. In these cases canard cycles are not allowed. 

o Case li- d < e, Case 2\. d = e and Case 3i. d > e: The points of 
S outside the interval {d,e) belong to Si. The points inside this interval, 
when it is not degenerated, belong to S3 in Case li and to S2 in Case 3i. 
In both cases H{z) > for all z S (d, e). See Figure [T3l 



A<0 




) 


^ 





A>0 


% 




21^-^ 



















Figure 13. Cases li, 2i and 3i. 

o Case 4i. d < S, Case 5i. d = S and Case 61. d > S: The points of 
S outside the interval {d,S) belong to Si. The points inside this interval, 
when it is not degenerated, belong to S3 in Case 4i and to S2 in Case 61 . 
In both cases H{z) > for all z G {d, S). See Figure [Til 



A<0 




5i 



A = 

U 



^ 



A > 




Figure 14. Cases 4i, 5i and 61. 



In Cases 7i — 17i we assume that Y presents the behavior Y~^. 



FOLD-SADDLE BIFURCATION IN NON-SMOOTH SYSTEMS 



11 



o Case 7i. A < -P, Case 81. A = -/3, Case 9i. -/3 < A < -li/2, Case 
lOi. A = — /3/2 and Case lli. — /3/2 < A < 0; The points of S outside the 
interval (d, i) belong to Si. The points inside this interval belong to S3. 
The direction function H assumes positive values in a neighborhood of d, 
negative values in a neighborhood of i and ff(A/3/(l + /?)) = 0. So, by dll), 
the S— attractor P = (A/3/(l +/3),0), nearby (0,0), is the unique pseudo 
equilibrium. In these cases canard cycles are not allowed. See Figure [TSl 



A< 




/3/2 < A < 



Figure 15. Cases 7i - 111. 

o Case 12i. A = 0.' Since a = —1 and d = i,\i is straightforward to show 
that each point Q G (/i, z) belongs to a closed curve composed by an arc of 
X and an arc of y. So d = i is a S— center. See Figure [T6l 




Figure 16. Case 12i. 



o Case 13i. < A < /3/2, Case 14i. A = /3/2, Case 15i. /3/2 < A < /3, 
Case I61. A = /3 and Case 17i. A > j3: The points of S outside the 
interval (i, d) belong to Si and the points inside this interval belong to S2. 
The direction function H assumes positive values in a neighborhood of d, 
negative values in a neighborhood of i and H[\(i/{1 + /3)) = 0. So, by 
([6]), the S— repeller P = (A/3/(l + /3),0), nearby (0,0), is the unique pseudo 
equilibrium. In these cases canard cycles are not allowed. See Figure [T7l 

The bifurcation diagram is illustrated in Figure [TBI D 
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< A < /3/2 A = ^/2 /3/2 < A < /3 A = /3 A > /3 






Figure 17. Cases 13i - 17i 




Figure 18. Bifurcation Diagram of Theorem 1. 



4. Proof of Theorem 2 

Proof of Theorem 2. In Cases I2, 22 and 82 we assume that Y presents the 
behavior Y~ . In Cases 42, 82 and 62 we assume that Y presents the behavior 
Y^. In Cases 72 — 192 we assume that Y presents the behavior Y~^. 

o Case I2. d < e, Case 22. d = e, Case 32- d > e, Case 42. d < S, Case 
52- d = S and Case 62- d > S: Analogous to Cases li, 2i, 3i, 4i, 5i and 
61. 

o Case 72. A < -(3, Case 82. A = -/?, Case 92- -/3 < A < -/3/(l - a), 
Case IO2. A = — /3/(l — a) and Case II2. — /?/(l — a) < A < 0.' Analogous 
to Cases 7i — lli changing — /3/2 by — /3/(l — a) = —dist{h,i)/2, where 
dist{h,i) is the distance between h and i. The unique pseudo equilibrium 
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occurs in P = {p ,0) where 

(7) 

p-= ^i \M l-a){l + P) + \{l + a)+ 

2(a + 1) 

- V((l - a)(l + /3) + A(l + a)Y - 4/3(1 + a)(l + a + A(l - q))). 

o Case 122- A = 0; The points of S outside the interval {d,i) belong to 
Si and the points inside this interval belong to S3. The direction function 
H assumes positive values in a neighborhood of d, negative values in a 
neighborhood of i and H{pq,Q) = where p^ is given by ([7]) with A = 0. 
So P = (P(|,0) is a S— attractor. Since e = 0, it is easy to see that there 
is an arc ^i- of X connecting the points h and j. It generates a S— graph 
r = 7]''- U (72 U S U o"! of kind I. Since — 1 < a < — I+Eq, it is straight forward 
to show that the First Return Map r] = ify o i^x , where 



^x 



ipY 



s 

z = (x,0) 

(^,j)cs 

z = {x,0) 



(-x + 2A,0) 



and 



{h,i) CS 
x{ii + /3) - 2il 



/3-n 

has derivative bigger than 1 in the interval {h,d). By consequence, T is a 
repeller for the trajectories inside it and in this case canard cycles are not 
allowed. See Figure [T9l 



A = 0<A<n X = h 




Figure 19. Cases 122, 132 and 182. 



o Case 132. < A < ii.- The distribution of the connected components 
of S and the behavior of H are the same of Case 122 with P = {p~^,0) 
where p^ is given by ([7]). Since < A < ii, there is an arc ^y^ of X 
connecting the point j to a point fci G {h,d). Also there is an arc 7^ of 
Y connecting the point ki to a point /i € («,j). Repeating this argument, 
we can find an increasing sequence (fei)igN- We can prove that there is an 
interval / C {ki,d) such that r/' = (c/9y o cpx)' < 1- As P is a S— attractor, 
there is an interval J C {ki,d) such that t]' > 1. Moreover, there exists an 
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unique point Q G {k\,d) given hy Q = {{—if + X{ii + /3))//3,0) such that 
■q' = 1. By Q passes a repeller canard cycle T of kind I. See Figure [T9l 

o Case 142- X = ii: Every point of S belongs to Si except the point 
d = i. As in the previous case, we can construct sequences {ki)i^m and 
(/j)jgN- Since e = ii, we have that ki ^ d and U -^ d. So d is a non generic 
tangential singularity of repeller kind. In this case canard cycles are not 
allowed. See Figure [T9l 

o Case 152- ii < X < a/3/{l — a), Case I62. A = a/3/(l — a), Case 172- 
a/3/(l — a) < X < {3, Case I82. X = (3 and Case 192. X > f3: Analogous to 
Cases 13i — 17i changing j3/2 by a/3/(l — a) = —dist{i,j)/2. The unique 
pseudo equilibrium occurs in P = (p~,0) where p~ is given by ([7|). 




Figure 20. Bifurcation Diagram of Theorem 2. 



The bifurcation diagram is illustrated in Figure [20l 



D 



5. Proof of Theorem 3 

Proof of Theorem 3. In Cases I3, 23 and 83 we assume that Y presents the 
behavior Y~ . In Cases 43, 83 and 63 we assume that Y presents the behavior 
Y^. In Cases 73 — I93 we assume that Y presents the behavior Y~^. 

o Case I3. d < e, Case 2^. d = e, Case 83. d > e, Case A3, d < S, Case 
53. d = S and Case 63. (i > S: Analogous to Cases li, 2i, 3i, 4i, 5i and 
61. 

o Case 73. A < -(3, Case 83. A = -/?, Case 93. -(3 < X < -(3/{l - a), 
Case IO3. A = — /3/(l — a) and Case II3. —13/(1 — a) < X < ii: Analogous 
to Cases 7i — lli changing — /3/2 by — /3/(l — a) = —dist{h,i)/2. The unique 
pseudo equilibrium occurs in P = {p~ , 0) where p~ is given by ^ . 

o Case 123. A = ii: Analogous to Case 142 except that here d is an 
attractor, i.e., there is a change of stability. See Figure [211 
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Figure 21. Cases 123, 183 and I33. 



o Case 183. ii < A < 0; Analogous to Case 182 except that there is a 
change of stabihty on P = (p~,0), which is a S— repeller, and on F, which 
is an attractor canard cycle of kind I. See Figure ETJ 

o Case 143. A = 0; Analogous to Case 122 except that occurs a change 
of stability on P = (p~,0), which is a S— repeller, and on F, which is an 
attractor for the trajectories inside it. See Figure [211 

o Case 153. < A < a/3/(l - a), Case I63. A = a/3/(l - a), Case ITs. 
a/3/(l — a) < X < j3, Case I83. A = /3 and Case 192. A > /3.- Analogous to 
Cases I81 — 17i changing /3/2 by a/3/(l — a) = —dist{i,j)/2. The unique 
pseudo equilibrium occurs in P = (p~,0). 




Figure 22. Bifurcation Diagram of Theorem 3. 



The bifurcation diagram is illustrated in Figure 



D 
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6. Proof of Theorem A 

Proof of Theorem A. Since in Equation ([5]) we can take a in the interval 
(— oo, 0), from Theorems 1, 2 and 3 we derive that this equation, with r = i, 
unfolds generically the (Invisible) Fold— Saddle singularity. 

Observe that the bifurcation diagram contain all the typical configurations 
and all the distinct topological behavior described in Theorems 1, 2 and 3. 
So, the number of typical configurations is 55 and the number of distinct 
topological behaviors is 11. Moreover, each topological behavior can be 
represented respectively by the Cases li, 4i, 7i, 12i, 13i, 122, 132, 142, I23, 
133 and 143. 

The full behavior of the three— parameter family of non— smooth vector 
fields presenting the normal form ([5]), with t = i, is illustrated in Figure 
[23] where we consider a sphere around the point {\,fj,,(3) = (0,0,0) with a 
small ray and so we make a stereographic projection defined on the entire 
sphere, except the south pole. Still in relation with this figure, the numbers 
pictured correspond to the occurrence of the cases described in the previous 
theorems. As expected, the cases 5i and 52 are not represented in this 
figure because they are, respectively, the center and the south pole of the 
sphere. D 



7. Proof of Theorem 4 

Proof of Theorem 4- Since X has a unique S— fold point which is visible we 
conclude that canard cycles are not allowed. 

In Cases I4, 24 and 34 we assume that Y presents the behavior Y~. In 
Cases 44, 54 and 64 we assume that Y presents the behavior Y^. In these 
cases, when it is well defined, the direction function H assumes positive 
values. 

o Case I4. d < e: The points of S inside the interval (d, e) belong to Si. 
The points on the left of d belong to S3 and the points on the right of e 
belong to S2. See Figure [Ml 

o Case 24. d = e: Here Si = 0. The vector fields X and Y are linearly 
dependent on d = e which is a tangential singularity. Moreover, it is an at- 
tractor for the trajectories of Z crossing S3 and a repeller for the trajectories 
of Z crossing S2. See Figure [Ml 

o Case 34. d > e: The points of S inside the interval (e, d) belong to Si. 
The points on the left of e belong to S3 and the points on the right of d 
belong to S2. See Figure [Ml 

o Case 44. d < S: The points of S inside the interval {d,S) belong to 
Si. The points on the left of d belong to S3 and the points on the right of 
S belong to S2. See Figure [25l 

o Case 64. d = S: Here Si = and S is an attractor for the trajectories 
of Z crossing S3 and it is a repeller for the trajectories of Z crossing S2. 
See Figure [25l 
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/3<0< \ >/3>0 

/3 = 

Figure 23. Bifurcation diagram of the (Invisible) Fold— Saddle singularity. 



A<0 






A = A >0 



^ 



34 



Figure 24. Cases I4, 24 and 84. 



o Case 64. d > 5".' The points of S inside the interval {d, S) belong to 
Si. The points on the left of S belong to S3 and the points on the right of 
d belong to S2. See Figure [25l 
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A<0 



44 



))^ • ■ ^ - -A— » 




54 



A = A >0 




^ 



64 



^-^i» - ^ ■ -^» 




Figure 25. Cases 44, 54 and 64. 



In Cases 74 — 184 we assume that Y presents the behavior Y~^ . 

o Case 74. d < h, Case 84. d = h and Case 94. h < d < i: The points of 
T, inside the interval {d, i) belong to Si. The points on the left of d belong 
to S3 and the points on the right of i belong to S2. The direction function 
H assumes positive values on S3 and negative values in a neighborhood of 
i. Moreover, iJ(/3A/(-l + /3)) = and the S-repeller P = (^A/(-l + /3),0) 
is the unique pseudo equilibrium. See Figure [26l 

o Case IO4. d = i: Here Si = 0. The vector fields X and Y are linearly 
dependent on the tangential singularity d = i. A straightforward calculation 
shows that H{z) = (1 — /3)/2 ^ for all z € S/{d}. So d = i is an attractor 
for the trajectories of Z crossing S3 and a repeller for the trajectories of Z 
crossing S2. Moreover, A = {(i}u47'U<T2U{5}UcriU/i(i is a S— graph of kind 
III in such a way that each Q in its interior belongs to another S— graph of 
kind III passing through d. See Figure | 




Figure 26. Cases 74 - IO4. 



o Case II4. i < d < j, Case I24. d = j and Case 184. 
of S inside the interval (i, d) belong to Si. The points on 
to S3 and the points on the right of d belong to S2. The 
H assumes positive values on S2 and negative values in a 
Moreover, H{(3X/{-l + /3)) = and the S-attractor P = 
is the unique pseudo equilibrium. See Figure [271 

The bifurcation diagram is illustrated in Figure [28l 
behavior can be represented respectively by Cases I4, 24, 
II4. 



j < d: The points 
the left of i belong 
direction function 
neighborhood of i. 

= (/3A/(-l + /3),0) 

Each topological 
44, 54, 74, IO4 and 

D 
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0< A</3 




Figure 27. Cases II4 - 184. 




Figure 28. Bifurcation Diagram of Theorems 4, 5 and 6. 



8. Proof of Theorem 5 

Proof of Theorem 5. The direction function H has a root Q = (q, 0) where 

(8) 

■((-l + a)(l-/3)-A(l + a)+ 



2{a + l 



+v/((-l + a)(l -13)- A(l + a))2 + 4/3(1 + a)(l + a + A(-l + a))). 

Moreover, i:f assumes positive values on the right of Q and negative values 
on the left of Q. Note that when a^-— IsoQ^— 00 under the line 
{y = 0} ^iicl it occurs the configurations showed in Theorem 4. 

In Cases I5, 25 and 85 we assume that Y presents the behavior Y~. In 
Cases 45, 65 and 65 we assume that Y presents the behavior Y^. In Cases 
75 — 185 we assume that Y presents the behavior Y^. 

o Case I5. d < e, Case 2^. d = e, Case 85. d > e, Case 4^. d < S, Case 
55. d = S and Case 65. d > S: Analogous to Cases I4, 24, 84, 44, 84 and 
64 respectively, except that here it appears the S— saddle Q on the left of d 
and e or S. See Figure [29l 
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A < (1 + a)/3/(l - a) A = (1 + a)/3/(l - a) A > (1 + a)/3/(l - a) 





Figure 29. Cases I5, 25 and 85. 



o Case 75. d < h, Case 85. d = h, Case 95. h < d < i: Analogous to 
Cases 74 — 94, except that here it appears the S— saddle Q on the left of d 
and i. Here P = (p, 0) where 
(9) 



P 



2{a + i 



-((-l + a)(l-/3)-A(l + a)+ 



V((-l + a)(l - /3) - A(l + q))2 + 4/3(1 + a)(l + a + A(-l + q))). 



o Case IO5. d = i: Analogous to Case IO4, except that here appear the 
S— saddle Q on the left oi d = i. 

o Case II5. i < d < j, Case I25. d = j and Case 185. j < d: Analogous 
to Cases II4 — 184, except that here it appears the S— saddle Q on the left 
of d and i. 

The bifurcation diagram is illustrated in Figure [251 Each topological 
behavior can be represented respectively by Cases I5, 25, 45, 55, 75, IO5 and 
II5. □ 



9. Proof of Theorem 6 

Proof of Theorem 6. The direction function H has a root Q = {q, 0) where 
q is given by ([8]). Moreover, H assumes positive values on the left of Q and 
negative values on the right of Q. Note that when a^— IsoQ^-oo under 
the line {y = 0} and it occurs the configurations showed in Theorem 4. 

o Case 1q. d < e, Case 2q. d = e, Case 3q. d > e, Case Aq. d < S, Case 
5q. d = S and Case 6q. d > S, Case 7q. d < h, Case Be- d = h, Case 9q. 
h < d < i, Case lOe- d = i, Case lie- i < d < j, Case 126- d = j and Case 
ISe- j < d: Analogous to Cases I5, 25, 85, 45, 65, 65, 75, 85, 95, IO5, II5, 
125 and I85 respectively, except that here the S— saddle Q takes place on 
the right of d, e, S and i when these points appear. 

The bifurcation diagram is illustrated in Figure [28l Each topological 
behavior can be represented respectively by Cases Ig, 20, 46, Se, 76, 106 and 
llfi. D 



FOLD-SADDLE BIFURCATION IN NON-SMOOTH SYSTEMS 



21 



10. Proof of Theorem B 

Proof of Theorem B. Since in Equation ([5]) we can take a in the interval 
(—1 — eo, —1 + eo) we conclude that Theorems 4, 5 and 6 prove that this 
equation, with r = u, unfolds generically the (Visible) Fold— Saddle sin- 
gularity. Its bifurcation diagram contains all typical configurations and all 
distinct topological behavior described in Theorems 4, 5 and 6. So, the 
number of typical configurations is 39 and the number of distinct topologi- 
cal behavior is 21. Moreover, each topological behavior can be represented 
respectively by the Cases I4, I5, le, 24, 25, 2q, 4^, 4^, Aq, 54, 85, 5q, 7^, 7^, 
76, IO4, IO5, lOe, II4, II5 and lie- 




A = 



A<0 



/3<0 < \ >/3>0 

/3 = 

Figure 30. Bifurcation diagram of the (Visible) Fold— Saddle singularity. 



The full behavior of the three— parameter family of non— smooth vector 
fields presenting the normal form ([5]), with t = v, is illustrated in Figure 
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[301 where we consider a sphere around the pomt (A, /z,/3) = (0,0,0) with a 
smah ray and so we make a stereographic projection defined on the entire 
sphere, except the south pole. Stih in relation with this figure, the numbers 
pictured correspond to the occurrence of the cases described in the previous 
theorems. As expected, the cases 64 and 65 are not represented in this 
figure because they are, respectively, the center and the south pole of the 
sphere. D 
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